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1. INTRODUCTION 

Batteries are used in almost every part of our lives. Increasingly, it is finding its utility in new 
application areas such as hybrid electric cars, personal mobile devices, micro energy harvesting system and 
so on. There are extensive research efforts in many areas of the battery technology [1]-[5]. One of the 
important area is in the modeling of the battery because this allows researchers, designers and engineers to 
evaluate various design choices in a timely and cost effective manner. However, most of the work has been 
focused on Li-ion batteries due to their application in hybrid electric cars [6]-[11]. On the other hand, in the 
rapidly emerging areas of wireless sensor network (WSN) and internet of things (IoT), there has been less 
attention for the smaller capacity nickel metal hydride (NiMH) batteries. Therefore, this work endeavours to 
address this gap by developing a simple but accurate NiMH battery model. 

WSN and IoT are experiencing an explosive growth due to their wide applicability across many 
industries. In order to enable the continuous and perpetual operation of WSN and IoT, micro energy 
harvesting system has emerged as a promising solution. A micro energy harvesting system consists of energy 
harvesting, energy storage and power management subsystems. The energy storage subsystem usually 
comprises a battery as a store of energy as renewable energy source. Since this energy source supplies in a 
non-constant manner, the system operation can be disrupted. Therefore, the battery is a key component in the 
micro energy harvesting system. 

The ability of a battery model in simulating its behaviour under various conditions not only reduces 
time but also cost when compared to building a hardware prototype. Various tests can be performed and their 
results obtained in a timely manner [12]. In the literature, a large variety of models have been proposed, 
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which include electrical, electrochemical, analytical, and stochastic models. Among them, the electrical 
models are preferred due to its good trade-off between accuracy and simpler parameterization effort. 
Furthermore, these models also take into account both the complexity and non-linearity of the battery [13]. In 
contrast, electrochemical models require the solving of partial differential equations, which is much more 
complex. Therefore, the electrical models offer good usability and accuracy [14]. 

The electrical battery model comprises of an equivalent circuit with voltage source, resistance and 
capacitance. Its variants include one resistive (1R), one resistive-capacitance model (1RC) and two resistive- 
capacitance (2RC) models. The 1R model is represented by only series resistance with a voltage source. On 
the other hand, IRC and 2RC have an additional single RC branch and double RC branch, respectively. 
Among them, the 1RC offers the best trade-off between accuracy and parameterization effort [15]. 

The accuracy of the model is highly dependent of the parameterization, which is basically to 
determine the values of the resistors and capacitors that should be used in the corresponding model. Two 
main techniques used for parameterization are the pulse measurement and optimization. A rectangular current 
pulse is applied to obtain the voltage response of the battery in both techniques. For the optimization 
technique, the values of resistance and capacitance are optimized by matching the produced voltage response 
against the measured voltage response. This process has a higher complexity and involves significant efforts. 
In the pulse measurement technique, the values of the resistance and capacitance are obtained through 
mathematical equations from the measurement of the voltage response [16], which is much simpler. Among 
the pulse measurement based parameter extraction techniques, one of the best was proposed Einhorn et al. in 
[15]. It consists of simple equations that can be applied with the pulse measurement procedures. The 
parameter extraction technique has been applied and tested on Li-ion battery, with very high accuracy. 
However, it has not been applied on the NiMH battery. 

Although limited, there are several recent works in NiMH battery modeling. Cruz-Manzo et al. and 
Zhu et al. in [17], [18], electrochemical models were proposed with the objectives of obtaining new 
understanding in the interpretation of battery electrochemical mechanisms and to study the battery capacity 
degradation effects. These models are complicated and its integration with other electronic components’ 
models is also not straightforward. An investigation into the tradeoff between accuracy and simplicity was 
undertaken in Fotouhi et al. [19] based on the electrical model and optimization-based parameterization. The 
study yielded very accurate results in terms of state of charge (SoC) and state of health (SoH). Meng et al. in 
[20], an accurate SoC estimator for design of battery management system (BMS) was proposed for battery 
lifetime extension. It worked very well in the specific application of electric boats that were to be used for 
scenic tours in caves. A enhancement to an existing analytical based battery model by taking into account the 
temperature effect was undertaken Rodrigues et al. in [21]. It provided an accurate estimation of the battery 
lifetime at different temperatures. However, in these works, the detailed voltage response results were not 
provided. 

Therefore, this paper will present the IRC NiMH battery model with the simple parameterization 
technique. This is followed by its implementation in a physical and dynamic simulator. Experiments with 
NiMH batteries using both current pulses and real wireless sensor node current profiles will be performed to 
test and validate the model's accuracy. The percentage difference between the model and experiment is less 
than 2%. The proposed parameterization technique shows a high accuracy with NiMH battery compared to 
the previous works, which only tested with Li-ion batteries. Following is the organization of the paper. The 
methodology for IRC battery model parameterization is presented in section 2. The results are shown and 
discussed in section 3, while section 4 concludes the paper. 


2. METHODOLOGY 
2.1. 1RC electrical battery model 

The battery modeling along with its parameter extraction technique is presented. The model is 
relatively simple and requires minimum parameters from the datasheet. IRC electrical battery model consists 
of a voltage source V, series resistance R, and a parallel RC branch connected in series is shown in Figure 1. 
The model is capable of predicting the battery performance and SoC estimation better than the typical 
analytical and electrochemical battery models. The model components are: 

— Vis the voltage source and it represents the open circuit voltage (OCV) of the battery 
— Ry is the series resistance of electrolyte. 
R», Cp is the RC network and represents the transient response of the battery electrodes. 

Parameterization is to extract the values of model parameters such as Rs, Rp and Cp. The IRC battery 
model is dynamic in nature and the parameter values are varying with the change in SoC. Therefore, the 
important part is to get the correct parameters of the model when they are varying with SoC. The IRC 
electrical battery model parameters are: 
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— Cy, the nominal capacity of a battery. 

— Rs, the series resistance of electrolyte at SoC=0. 
— Rp, the resistance of electrodes at SoC=0. 

— Cp, the capacitance of electrodes at SoC=0. 





Figure 1. 1RC battery model 


The SoC can be obtained by integrating the current over time. The SoC varies between O and 1. The 
SoC vs OCV curve provides the change in SoC of the battery with voltage. It is determined when there is no 
load attached to the battery and is shown in Figure 2. To extract the parameters of IRC battery model, an 
experiment procedure to obtain the pulse measurement is required. Firstly, the battery is fully discharged to 
0% SoC followed by a rest period. The rest period is to ensure the accurate measurement of the OCV after 
obtaining the thermodynamic stability [22]. After the rest period, a current pulse is applied for a certain 
duration to obtain the voltage response of the battery. The cycle of charging and resting period is repeated 
until the SoC reaches 100%. This whole procedure will provide the SoC vs OCV curve. Alternatively, this 
procedure can also be performed on a fully charged battery by applying discharging current pulses with the 
following rest periods. Hence, after obtaining the SoC vs OCV curve, the parameters can be extracted using 
the equations provided Hentunen et al. in [23]. 
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Figure 2. SoC vs. OCV curve of NiMH battery 


2.2. Parameterization of the 1RC electrical battery model 

The manufacturer provides the nominal capacity, Cy, of the battery in the datasheet. The values of 
Rs, Rp and C, are obtained by applying a current pulse and measuring the voltage response. The 
parameterization technique is following the approach of [15]. The state of charge of the battery is represented 
in (1): 


SOC = CN- (1) 


CN 


where, Q is the charge of the battery at time to given as shown in: 
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t_0 
Q = J Is()dt (2) 
The series resistance of the electrolyte can be approximated with (3): 


R = 7) (3) 


where, / is the height of the pulse current and V; is the initial voltage response. While resistance of electrodes 
are shown in (4): 


_Yf 
Rp =- i (4) 


where, Vpis the final voltage response. Both V; and Vy are shown in Figure 3. The curve is obtained by 
applying a rectangular current pulse of 1.5A to the battery. 
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Figure 3. Battery voltage response 


The step current of the pulse measurement technique includes a series of charge or discharge cycles 
followed by the rest period. The width of the pulse is normally from 1 minute to 5 minute. In this paper, the 
l-minute pulse width is chosen to capture any abrupt change in the behaviour of the battery. A 1500 mAh 
NiMH battery is chosen to perform the tests. The step current of the pulse measurement technique is applied 
to the battery and the voltage response is obtained. The battery is initially discharged to 0% SoC and then 
subject to partial charge phase cycles until the SoC reaches 100%. After every charge cycle, a rest is given to 
the battery to stabilize the voltage. Therefore, at the end of one-hour cycle of rest period, the battery voltage 
is found stable to be considered a good estimate of the OCV. The SoC is determined based on the current 
gained or drawn from the battery at each cycle. This technique is also called coulomb counting. The 
capacitance of electrodes can be approximated as (5): 


a) 
Cp = Te (5) 


where, t is the time constant of the exponential curve and is determined by (6): 
T= Q(1-—1/e) (6) 


where e is the Euler’s number=2.71828. After obtaining the SoC vs. OCV curve by applying the step current 
of the pulse measurement technique, the parameters are extracted using (3)-(5). 


2.3. Implementation of 1RC electrical battery model 

The one diode PV panel model can be implemented through simple MATLAB equations. However, 
since the battery is dynamic, its model implementation requires both physical and dynamic simulators such as 
Simscape and Simulink. The model utilizes Simscape to build custom circuit elements such as source, 
resistor and capacitor which are dependent on SoC [24]. 
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The IRC NiMH equivalent circuit model system architecture is shown in Figure 4. The major 
building blocks are the signal builder, IRC equivalent circuit model, ideal temperature source, convective 
heat transfer, voltage sensor, and the scope. The signal builder provides the pulse charging current. The ideal 
temperature source basically represents an ideal source of thermal energy that is able to maintain the 
specified temperature of the system. The convective heat transfer block simulates the heat transfer by 
convection. Furthermore, the voltage sensor block converts the voltage measured into a physical signal that 
can be captured and displayed. 
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Figure 4. IRC equivalent circuit model system architecture in Simulink and Simscape 


The IRC model is required to generate the simulation results which are later compared with the 
experiment to verify the accuracy of the model. The 1RC equivalent circuit model has four variable circuit 
elements that are SoC dependent. These variable circuit elements are open circuit voltage Vəc, series 
resistance R, and RC network, R, and Cp. These variable circuit elements are implemented through lookup 
table. An example of the lookup table is shown in Table 1. The built-in Simulink thermal model is used to 
model battery temperature. It is assumed that the heating is primarily from internal resistance and the cooling 
is primarily via convection. The model can be used to simulate various operating conditions by just changing 
the temperature and current profiles. 


Table 1. Example of lookup table used in simscape model 


SoC OCV _ R,(Ohm) R, (Ohm) G,(P) 


0.2 1.24 2.13 1.54 6134 
0.4 L27 2.15 1:05 6207 
0.6 Loa 2.67 1.73 6398 
0.8 1.38 2.45 1.43 6254 
1.0 1.43 219 1.32 6243 


3. RESULTS AND DISCUSSION 

Validation of the model is carried out through experiments. The main component of the 
experimental setup is the ubiquitin like modifier activating enzyme 5 (UBAS) battery analyzer. The battery 
analyzer consists of charging/discharging and measurement subsystems. The former allows 
charging/discharging at constant voltage, constant current or constant power modes. The test cycle, which is 
user programmable, is used to control the charging and discharging of the battery. The latter includes a built- 
in measurement system that can allow simultaneous measurement of both the current and voltage. All of 
these subsystems are controlled by the battery analyzer and its associated software which is installed on a PC. 

In the experiments, the battery analyzer is used to apply the required current profile and measure the 
corresponding voltage response [25]. The NiMH battery can be charged or discharged with a completely 
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programmable current profile. The NiMH battery with capacity of 1500 mAh from Energizer is used in the 
experiment. The battery is tested to verify the model’s accuracy at room temperature. The block diagram of 
the experiment is shown in Figure 5. 

The step current pulses of 1-minute width are applied to both the NiMH battery and the model. By 
using (3)-(5), the parameter of IRC battery model can be obtained and apply to the Simulink. The 
comparison of the corresponding voltage responses obtained is shown in Figure 6 (a) and 6 (b). It can be seen 
that they are in excellent agreement. The difference between the two curves is less than 1.5%, which shows 
the high accuracy of the model. 
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Figure 6. Voltage responses obtained from experiment and model with, (a) 1 minute pulse, (b) 2 minute pulse 


Bulletin of Electr Eng & Inf, Vol. 10, No. 4, August 2021 : 1793 — 1802 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 1799 


The validation of the model is also performed for both 1 minute and 2-minute pulse width in order 
to test the accuracy of the model. The step current of 2 minute is applied to both NiMH battery and the model 
in order to validate its accuracy. The 2-minute pulse width is equivalent to 3.2% of the SoC of the battery. 
The percentage difference between the voltage response obtained from both model and NiMH battery is also 
less than 1.5% which is very accurate. 

Furthermore, the NiMH battery model is also tested with the real wireless sensor node current 
profiles obtained from [26]. The current profile of the WSN node includes the sleep, wakeup and ON 
conditions of both the micro-controller unit (MCU) and the radio. Both the current profiles during 
transmission and reception are considered. The current profile during transmission is shown in Figure 7 (a). It 
can be seen that the highest current consumption of the wireless sensor node occurs during the radio 
transmission mode. The current profile in reception is also similar to the transmission mode. There are 
sudden changes in the current profile of the wireless sensor node when transitioning from one state to 
another. The NiMH battery model should be able to allow for these sudden changes of the current profile and 
reflect them accurately in the output battery voltage. 

The current profiles are applied to the battery model and the corresponding output voltages of the 
battery are obtained. Furthermore, the same current profiles are also applied to the real NiMH battery. The 
comparison of the output battery voltages of transmission is shown in Figure 7 (b). The percentage difference 
between the output battery voltages of model and experiment is less than 2% which is highly accurate. This 
compares well with works on Li-ion batteries, where the accuracy reported ranges from 2%-5% [15], [22]. 
After testing and validating the model with the real current profile of the wireless sensor node, it can be 
concluded that the model can be adopted for any sophisticated current profiles of the wireless sensor 
networks (WSNs) and IoT devices. 
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Figure 7. These figures are; (a) current profile of wireless sensor node in the transmission mode, 
(b) comparison of voltage response obtained from NiMH battery and model 
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4. CONCLUSION 

IRC electrical battery model with simple parameterization for NiMH battery has been presented in 
this paper. The parameters are obtained through the pulse measurement technique which includes applying a 
series of charge cycles followed by the rest periods. Subsequently, the model is implemented in the physical 
and dynamic simulator. To determine the accuracy of the model, two step currents and the real current 
profiles of the wireless sensor node are applied to both real-life NiMH battery and the model. Compared to 
the previous works, which parameterized Li-ion batteries only, the resulting NiMH voltage responses 
obtained from the proposed model show that they are in very good agreement with difference of less than 
2%. It can be concluded that the presented model is simple and can be applied to any sophisticated real-life 
current profiles of the WSN and IoT, with excellent accuracy. Moreover, the validation of proposed model 
with different types of battery chemistry is also recommended. 
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